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We compute the correlation functions for antiferromagnetic and d-wave 
superconducting fermion bilinears in a generalized mean field type approxi- 
mation for the Hubbard model. For high temperature our explicit expressions 
show that homogeneous field configurations are preferred for these composite 
bosons. Below a critical temperature we find spontaneous symmetry breaking 
with homogeneous expectation values of the composite fields. Our results can 
be used to device a nonperturbative flow equation for the exploration of the 
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iow temperature regime. 
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1 Introduction 

In the attempt to come towards a better understanding of the properties and basic 
mechanisms of high temperature superconductors, the Hubbard model fl] has be- 
come one of the most studied models for systems with strongly correlated electrons. 
However, it has turned out that it is very difficult to extract any thermodynamic 
properties of this model in the interesting range of temperatures, couplings and car- 
rier densities even on a numerical level. What renders this problem so difficult is the 
fact that the most prominent physical degrees of freedom (such as antiferromagnetic 
or superconducting behavior) emerge as a complicated momentum dependence of the 
effective multi-fermion couplings induced by fluctuations and that different length 
scales are involved. We suggest that an easier understanding can be gained if the 
interesting degrees of freedom are included in a more explicit way. In an earlier paper 
we have presented a reformulation of the 2d Hubbard model, the so called colored 
Hubbard model f2|. In this reformulation, the interesting physical properties are de- 
scribed by bosonic degrees of freedom. These bosonic fields correspond to fermionic 
composite operators in the appropriate antiferromagnetic or superconducting chan- 
nels. In particular, this formulation allows us to extract the antiferromagnetic or 
superconducting properties from a calculation of the two point correlation functions 
of the bosonic fields instead of the more traditional investigation of quartic fermionic 
vertices. It thus becomes possible to apply standard calculation procedures to the 
problem. 

In this paper, we compute the one loop corrections to the bosonic propagators to 
get a first impression of how antiferromagnetic and superconducting behavior come 
into play when fluctuations are included. The intention to do this is twofold: First, 
we show that in the colored Hubbard model we are able to calculate a mean field 
type approximation for the correlation functions for composite operators analyti- 
cally. Our results are expected to be quantitatively reliable at high temperature if 
the interaction strength is not too large. Furthermore, they are expected to serve as 
a good qualitative guide even at relatively low temperature, where the complicated 
physics near the Fermi surface becomes dominant. We emphasize that the one loop 
calculation of the bosonic propagator accounts for contributions to the effective four 
fermion interaction which involve arbitrarily high powers of the coupling constant. 
Second, the calculations in this paper may serve as a starting point for a renormal- 
ization group analysis in the frame of the colored Hubbard model. Within an exact 
renormalization group approach || they allow us to motivate truncation schemes 
for the bosonic propagators by identifying the kinetic terms that emerge in our one 
loop calculation. 

Our work is organized as follows. In the next section, we present our results for 
the one loop correlation functions of the two most prominent composite degrees of 
freedom of the Hubbard model in the antiferromagnetic and d- wave-superconducting 



channels. We interprete the relevance of these results for the onset of spontaneous 
symmetry breaking. In the following two more technical sections, we give a brief 
review of the colored Hubbard model and present a mean field type approximation 
for the propagators of a whole set of bosonic fields. Additionally, in sect. 5 we give 
the high temperature limit, which has a particularly simple form and is well suited 
for the truncation ansatz for the renormalization group analysis mentioned above. 
Sect. 6 sketches briefly how nonperturbative flow equations can be derived from our 
results and sect. 7 contains our conclusions. 



2 Correlation functions for antiferromagnetic and 
superconducting behavior 

The Hubbard-model is defined for electrons on a lattice by the Hamiltonian 

ij,cr i 

where a ia and a^ are creation-/annihilation-operators for an electron at site i with 
spin a and obey the usual anticommutation relations {a ial dj :T } = 5ij5 aT . The prob- 
ability amplitude of an electron for tunneling from site i to site j is denoted by 
tij. We take a square lattice in 2 dimensions with ty = —t for neighboring lattice 
sites and otherwise. The interaction term ~ U mimics a screened Coulomb-like 
interaction, with m >ff = a i(T a ijff the number operator. 

Before turning to the technical details, we would like to illustrate our point by 
discussing our results for antiferromagnetic and superconducting properties of the Id 
Hubbard model on a square lattice. Obviously, antiferromagnetic or superconducting 
behavior is non-local in nature (e.g. antiferromagnetism emerges if spins of electrons 
situated on neighboring lattice sites are opposite in sign). We include this non- 
locality of the most interesting physical degrees of freedom by dividing up the lattice 
into square plaquettes of 4 lattice sites each, which we enumerate clockwise. In 
particular, we define electron-hole and electron-electron bilinears 



a(x) = ^*(2?)r0i(x) - il)l(x)Ti>2{x) + ^(x)t^3(x) - ii)l{x)fi> 4 (x) 
d(x) = z[$l(x)t 2 $ 2 (£) - ^2{x)t 2 ^{x) + i)z{x)r 2 ^i{x) - $ 4 (z)t 2 $i(x)], 



(2) 



where x = (2cm, 2am), n, m e Z, is the position vector of the corresponding plaque- 
tte (cf. fig. ([!])), with a the lattice distance. Nonzero expectation values of a or d 
correspond to antiferromagnetic and (d-wave)-superconducting states, respectively. 

In [[| we have shown that it is possible to rewrite the original partition function of the 
Hubbard model in a partially bosonized form by introducing a set of fermion bilinears 





Figure 2: The one loop correction Ar^ 2 ^ to the bosonic kinetic term in the effective 
action for the boson a and the boson d at high temperature, for T/U = 10, t/U = 1, 
h 2 /U = 10. 



similar to the two we wrote down explicitely in eq. (@). The full set of fermion 
bilinears is given in appendix ||. Additionally, to preserve translational invariance 
on the original lattice, we introduce a color index for the bilinears d c (x), d c (x), c = 
1 ... 4. Different values of c indicate a shift of the corresponding bilinear in position 
space by one lattice unit a. The exact definition is provided in appendix p. In the 
partially bosonized reformulation of the Hubbard model the bosonic fields a c , d c , . . . 



couple to their fermion-bilinear counterparts a c , d c 



and therefore describe the 



intended degrees of freedom explicitly. Before fluctuations are included the effective 
bosonic interactions are purely quadratic. The inverse bosonic propagators contain 
only the bare mass term which is the same for all bosons. It is of great interest to 
investigate the behavior of the quadratic (kinetic) bosonic terms once fluctuations 
are included. They describe the propagation of the corresponding composite degrees 
of freedom. Possible instabilities and the onset of spontaneous symmetry breaking 
are indicated by zeros of the quadratic terms. 

In this paper we have calculated the one loop corrections to each of the bosonic 
propagators in our partially bosonized theory. They correspond to a "mean field 
type" analysis where the effects of fermion fluctuations are included whereas boson 
fluctuations are omitted. All our results are for half filling (// = 0). In this section we 
concentrate on appropriate averages over the position index c (see eq. (|25|) ) which 
we denote by a and d. The results can be expressed in terms of the quadratic piece 
in the effective action T for bosons which reads in momentum space (we take the 
lattice spacing a = 1/2) 
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d* (m, k) ( 4tt 2 + Ar f (m, k) ) d(m, k) 





Figure 3: The one loop correction AT^ to the bosonic kinetic term in the effective 
action for the boson a and the boson d at low temperature, for T/U = 0.1, t/U = 1, 
h 2 /U = 10. We observe large negative values for k = (0, 0) indicating instabilities. 



+-a(-m, -k) Uir 2 + AT {2 \m, kj) a{m, k)\ + • ■ 



(3) 



For convenience, we suppress the dependence on the Matsubara frequency m in the 
following. 

The inverse bosonic propagators T^' a (k) = 

An 2 + Ar^(fc) are directly related to mo- 
mentum dependent four fermion interac- 
tions of the form 




(2), 



d*(k) (T^(k) 

a(-k) (ri 2 >(£) 



' d(k), 

-l _, _, 

a(k), 



(4) 



which arise from the exchange of the com- 
posite boson (see eq. (|T9"D). Here d(k), d(k) 
are the Fourier transforms of color aver- 
ages of the bosons in eq. (0) (see eq. (p5|)). 



Figure 1: The introduction of plaque- 
ttes and the enumeration of the colors. 



(2) 

In particular, the vanishing of T B for a 



given momentum k corresponds to a diver- 
gence of the four fermion interaction in the 
corresponding channel. In the bosonic lan- 
guage this instability is easily interpreted: 
In addition to the quadratic terms (0) the fluctuations also induce higher order 
bosonic interactions. They stabilize the effective scalar potential for large values of 
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Figure 4: The inverse propagator at zero momentum T a (0) as a function of tem- 
perature. The uncertainty of our mean field approximation is demonstrated by two 
different effective couplings h a . The left plot is for h 2 a /U = 10, the right plot for 
h 2 a /U = 40. We have indicated the lower bound for the critical temperature by a 
circle. 

a and d (see 0). A negative quadratic term therefore implies that the minimum 
of T occurs for nonvanishing values of the corresponding bosonic field and reflects 
spontaneous symmetry breaking. 

(2) 

In thermal equilibrium the inverse propagators T B depend on the temperature T. 
They typically decrease for decreasing T. The critical temperature T c for a phase 
transition to an antiferromagnet or a superconductor is equal to (for second order 
transitions) or above (for first order transitions) the temperature for which T B (T) 

(2) 

becomes negative first. The field for which T B becomes negative first is likely to 
indicate the order parameter for the low temperature phase. Similarly, the value of 
the momentum k for which T B reaches zero first tells us about the preferred spatial 
dependence of the order parameter. If in the vicinity of the phase transition the 
minimum of T B (k) occurs for fc = 0we expect a homogeneous bosonic expectation 
value. Such a scenario would lead to a considerable simplification since the dominant 
nonlocality of our problem would be absorbed in the definition of appropriate bosons 
as nonlocal fermion bilinears (^). 

Our results for the inverse bosonic propagators 1^ , {B = a,d), are shown in figs 
||-§]. In figs. H] and ^| we have plotted T B as a function of the bosonic momentum 
k for different temperatures. We choose the Matsubara frequency u m = of the 
bosons which is the dominant term in the propagator. The parameter h 2 /U reflects 
a freedom in the choice of bosonization and will be explained in section |3|. 

What can we learn from these plots? First of all, observe that for high temperature 
(fig. |2|) we get a simple momentum dependence. In this limit the one loop corrections 
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are of the form ~ cos 2 ^ cos 2 ^ for a and ~ (cos 2 ^ + cos 2 ^) for <i, namely 

Q/,2 

Ari 2 )(A;) = --^[cos(fc 1 /4)cos(fc 2 /4)] 2 , 

S (5) 

Arf(fc) = -^[cos 2 (^/4) + cos 2 (A; 2 /4)]. 

This has the nice consequence that the loop corrections in the high temperature 
limit do not change sign when varying parameters or momenta^]. The periodicity 
4tt in hi reflects the original lattice on which the Hubbard model is formulated. 
The momentum range of the bosons, however, is given by — n < ki < 7r so that 
the propagator for each boson has only one single minimum at k — 0. When T is 
lowered, the physics near the Fermi surface comes into play (in our formulation, 
the Fermi surfaces are given by k\ = ±fc2 m od27r) and the momentum dependence 
becomes more complicated (figM). 

We find that the minimum of the one loop corrections to the bosonic term in the 
effective action is situated at k = for all values of the temperature. Remarkably, the 
bosonic propagator (47r 2 + Ar^ ) _1 diverges at low temperature for small values of 
\k\. This has important consequences once the bosonic fluctuations (omitted so far) 

(2) 

will be included. The propagation of momentum modes near zeros of T B is strongly 
facilitated. One therefore expects that the physical behavior of the system should 
be dominated by fluctuations around such zeroes. Our result for the momentum 
dependence of the bosonic propagator is quite pleasing, because momentum modes 
with k « correspond to (almost) homogeneous bosonic field configurations in 
position space. 

For the minimum of T B (k) at k = the relevant parameters for the onset of instabil- 
ity are the mass terms M B = T B (k = 0). Negative values of M B indicate instability 
of the "symmetric" state without antiferromagnetism and superconductivity. The 
one loop result for the mass term can also be inferred from 0: (q = 2tcos(qi/2)) 

Ml = ^-2hlf J^ E tanh(^ 1+£ e 2 )) . (6) 

J-n (2vr) 2 £ ^ 1 ci + ec 2 

In fig.|] we show the temperature dependence of the mass term M 2 = Ta (0). For our 
parameters the change of sign of M 2 indicates a second order phase transition. We 
have indicated the critical temperature with a circle in fig.|]. This result coincides 
with the diagrams in for p = (or \x = 0) which were derived by using a mean 
field approximation with constant homogeneous background fields. 



4 This feature will facilitate the regularization of the bosonic propagators in a renormalization 
group approach. 



In the same way as for a and d we have calculated the one loop corrections to each of 
the bosonic propagators in our partially bosonized theory [fj]. As we have included 
16 different boson species and each boson carries an additional color index, we get 
a symmetric matrix of 64(64 + l)/2 different kinetic terms in the effective action. 
It turns out that the inverse propagators become minimal for homogeneous field 
configurations. Thus it is tempting to analyze first the matrix entries in the limit of 
homogeneous fields. We find that under the assumption of homogeneous field con- 
figurations most of the entries vanish and nonvanishing mixings between different 
bosons occur mainly between bosons of the same species but different color (there are 
some nonvanishing contributions mixing complex bosons of different species which 
we will discuss below). Diagonalising the 4x4 propagator matrices in color space 
for a given boson species allows us to rewrite the 64 x 64-propagator matrix of the 
bosons in a block diagonalised form. For general non homogeneous field configura- 
tions we find that the propagator matrix remains diagonal in color space. However, 
additional non vanishing terms occur in the propagator matrix mixing bosons of 
different species and same color. For this set of bosons with propagator matrices 
which are diagonal in color space we have calculated the one loop correction also for 
the propagators at nonvanishing momenta which corresponds to the effective action 
for non homogeneous bosonic field configurations. We find that the momentum de- 
pendence of Ar# in the high temperature limit is given by some linear combination 
of cos 2 ^r, cos 2 ^p an d cos 2 ^ cos 2 ^ also for the other bosons. The discussion above 
thus generalizes to the other bosons as well. 

In conclusion, the explicit notion of nonlocal degrees of freedom in the colored Hub- 
bard model yields a very powerful tool to discuss spontaneous symmetry breaking 
in a very simple and illuminating way. Although our one loop results for strong 
interactions will not produce a quantitatively correct picture at low temperature, 
they nevertheless give a first impression of how the composite fields associated with 
antiferromagnetism and superconductivity propagate. 



3 The colored Hubbard model 

The colored Hubbard model is a reformulation of the Hubbard model for electrons 
on a square lattice with next neighbor hopping^]. Let us briefly review the basic 
relations in the colored Hubbard model as described in 0. The partition function^] 
reads 

Z = exp {- ^Jt } f D ^ D ^ m * mD{v ex P { " ( S kin + S Y + Sj) } (7) 



5 It can also be used for a generalization of the Hubbard model M 

6 See appendix [A] for details of the Fourier transforms. We set the lattice distance of the coarse 
lattice to unity: 2a = 1. 



with 

s kin = y £{r a (Q)Pi(Q)MQ) 



Q 



71" 1 

+7T 2 E u} c {Q)up c {Q) + y E ™-yc(-Q)w-/c(Q) ) 

7 

S Y = - J2 S(Q-Q'-Q") 

QQ'Q" 

+ E^(Q)€(QOC c W', -or)U-QT) 

7 

^ = E { " E IW)W) + ^c(Q)^ c (Q)] - E L lc {-Q)w lc {Q) 

Q /3 7 

-[^(Wa(Q)+%(Q)^(<5)]} (8) 

and inverse classical fermion propagators P^ 

P*(Q) = iu n -2t(cos( qi /2)^ °) + cos(g 2 /2) ( ° £ )) . (9) 

The momentum vector Q = (u n , q) involves the Matsubara frequencies u n = 2irnT 
with n integer for bosons and half integer for fermions. The action S = Skin + Sy + Sj 
describes a coupling of electrons and holes denoted by Grassmann variables ip, ip* 
to bosonic fields u, u* and w. The temperature T is a free parameter^. The original 
four fermion interaction of the Hubbard model is encoded in the vertex factors V 
which are proportional to the Yukawa couplings h with h 2 ~ U. It can be recovered 
by solving the field equations for the bosons as functionals of the fermions. The 
sums over j3 run over the set of complex bosons, the sum over 7 over the set of real 
bosons. These sets are listed in appendix [B] and defined more precisely in . We 
use a and b for fermion color and V are also matrices in spinor space. Respectively, 
c denotes bosonic color. The partition function eq. (0) reduces to the one of the 
Hubbard model if the Yukawa couplings are chosen appropriately (cf. appendix fB|). 
The explicit vertices corresponding to the Hubbard model eq. (P are shown in 
appendix y. 

The effective action r[^,<^] in the partially bosonized formulation is introduced by 
the usual Legendre transform 



T$,<p\ = - In Z[fj, K] + ffjj>+ I Kip (10) 



7 A nonvanishing chemical potential [i could be included in the source of p as L pc = L' + \i 



where fj and K are expressed in terms of ip and p by 



V> 



5\nZ 
Sfj 



p = 



5\nZ 
5K 



(11) 



Here we use a collective notation for fermionic fields ip and sources fj {ip includes ip, 
ip*, fj includes rj, rj*) and for bosonic fields tp and sources K (tp includes the real and 
complex bosons). The mapping of our results for the partially bosonized effective 
action to the generating functional r^ [ip] of the one particle irreducible correlation 
functions within the original fermionic theory is straightforward ||. In presence of 
arbitrary bosonic sources K we may define 



TW$,<p,K\ = T$M- j K<p. 



The field equations for p 



S V (K) 



5p 



(12) 



(13) 



i>,K 



is solved by tpo[ip,K]. The effective action in the purely fermionic language is ob- 
tained by insertion of the solution 



This is easily verified by the fermionic field equation 



(14) 



£r«o 



5ip 



S V [K) 



K 



Sip 



5p 



+ 



K S <P 



V 



(15) 



ip,K 



where the rhs has to be evaluated for tp = p so that the second term vanishes. 
In the approximation of this paper 



-IpAiP^ABtpB + 2^( r B )<tt<Pt ~ Vab^A^bPv 



(16) 



the classical solution of ( pT3|) for K — reads 

fOa = {T i ^)~ T VAB,ri ! A'4'B 

and we obtain the lPI-generating functional in the fermionic language 

(2) 

Our computation of T B therefore amounts to an approximation to the 1 PI- four 
fermion vertex which can be compared with the results from other methods H. 



(17) 



(18) 
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(2) 

Finally we would like to stress the simple connection between 1^ and the fermionic 
4-point Green functions. First note that T B = (B(x)B(y)} , where B E 

{a, d, . . .}. The right hand side is the inverse of the the connected bosonic 2-point 
function. From the explicit expressions for the original fermionic partition function 
and our partially bosonized partition function |2|, one finds 

(B A (X)B B (Y)) conn . = ^8 AB 8(X -Y) + j^~ 2 [(B* A {X)B B {Y)) - (B A (X))(B B (Y))) 



(19) 



where h A , h B are the Yukawa couplings between B and B. As a second order phase 
transition into the antiferromagnetic or d-wave superconducting phase occurs for 
diverging correlation length, we expect T B to vanish when such a phase transition 
takes place. 



4 One loop corrections 

In the partially bosonized formulation eq. (0) the action S = Skin + Sy + Sj is purely 
quadratic in the fermionic variables (in absence of the fermionic sources, rj = rj = 0). 
The fermionic part of the functional integration is therefore Gaussian and can be 
performed in a standard way. For this purpose it is convenient to define 

*<» -(£$)■ (20) 

so that the fermionic part of the action takes the form 

^4E*'W'^*")' P = Po-AP (21) 



Pn 



-(P*) T \ A p_(c 


-A T 
B 


A = V w (Q',Q")w(Q'-Q"), 




B = 2V u (Q',Q")u(Q' + Q"), 




C = 2V u *(Q',Q")u*(Q' + Q"). 





The interaction terms A, B, C are to be interpreted as matrices with momentum 
labels Q', Q", fermion colors a, b and spinor indices inherited from the vertices V. 
Summation over the boson species and their indices is implied in the products Vw 
etc. 

11 



Omitting the bosonic fluctuations, the effective action| becomes r = S + AT, where 
the arguments of S are now the "background fields" ip, u, w. The bosonic fields u, 
w appear in the contribution from the fermionic fluctuations Ar = — l/2TrlnP 
since A, B and C are linear in these fields. Therefore Ar is a functional of the 
bosonic fields. We emphasize that Ar accounts for the full one loop correction! to 
the bosonic part of V since S does not contain purely bosonic vertices. 

We next expand the loop correction Ar in the number of boson fields 



Ar = — TrmP = --Trln[P (l-P - 1 AP)] 



-~ (TrlnPo - Tr(P l AP) - ^(P^AP) 2 + • • ■) 



(22) 



The first field independent term is discarded and the second linear "tadpole" term 
does not contribute to the propagators of interest here. It vanishes for half filling 
(/i = 0). The bosonic propagator corrections are described by the third term. 

Ar 2 = ^[^ 7C (-g)tr{(p^)- 1 (Q')^(Q / ,Q + Q , )(^)" 1 (Q + Q') 



QQ' 

r W„J 



V/(Q + Q',Q')jw, /c iQ) 

-2n; c (g)tr{(p^)- 1 (g')^(Q^Q-Q / )((^)~ 1 ) T (Q-Q , ) 

vf{Q-Q' 1 Q'))u pic> {Q)\. (23) 



Here tr refers to a trace over color and spinor indices and V c denotes the matrix 
[Vc]ab = V a b, c - We note that Ylo i n ^1^ involves a summation over momenta and 
Matsubara frequencies m. We extract the inverse propagator matrices (K = (uo m , k)) 



= ^tr{(p^)' i (Q)^r(Q,^+Q)(^)^ 1 (^+Q)v: c y(^+Q,Q)} 

Q 
^lu 0lc XK) (24) 

= -2^tr{(P^)- 1 (Q)Kc' 3 (Q,^-Q)((P^)- 1 ) T (ir-Q)y c U ?'(K-Q,Q)} 



8 Field independent additive constants are neglected in T. 

9 One loop corrections to V which involve fermion fields are not described by Ar. 
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These one loop expressions 

K + Q 




K 

ir - - _, ^ 1| I ' I - - — ^ ■ - u 

K 



involve a momentum integration over (f and a sum over fermionic Matsubara frequen- 
cies n, with K = (u m , k), Q — (u n , q). The computation of the inverse propagators 
r( 2 >(K) is the aim of this note. We note that the one loop expression Ar( 2 ) is ~ h? 
and therefore ~ U. Retranslating our result to the fermionic language (eq. (£D) it 
corresponds to a resummation involving arbitrarily high powers of U. 



4.1 Representation with respect to translations 

When calculating the loop corrections for a given boson-species one observes that 
different colors are mixed. It is favorable to use color combinations that render the 
propagator matrix diagonal. This can be achieved by using combinations which are 
simple representations with respect to translations in x- and y- direction, i.e. the 
color-combinations 



B 1 = \{B X + B 2 + B 3 + S 4 ), B 2 = \{B X -B 2 + B 3 - B A ) 
B 3 = \{B X + B 2 ~B 3 - S 4 ), B 4 = \{B X -B 2 - B 3 + B A ) 



(25) 



To motivate this choice, note that for homogeneous bosonic fields (B c (x) = B c (x + 
e x / y ) the group of translations by a = 1/2 in the x- and y-direction is isomorphic 
to G = Z 2 x Z 2 . G is Abelian, has order 4 and therefore 4 irreducible (necessarily 
one dimensional) representations given by the parity of the two Z 2 factors. Four 
linear independent basis functions for the irreducible representations of G which 
yield eigenvalues ±1 are easily written down and are exactly the Bi. In particular, 
a homogeneous field B\ is invariant under translations by a. It turns out that the 
combinations eq. ( p5| ) diagonalize the propagator matrix in color space even in the 
case of nonvanishing momenta. This requires that the Fourier transforms are chosen 
as in appendix |A|. 

Writing eq. fl25D as B a = M ab B b /A 1 we define the vertices in the new basis as 
Vj = M ah V^. With M = M T , M 2 = 4, the Yukawa interaction S Y then looks the 
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same as in eq. ([/]), with B — > B and V B — > V B . Hence the above calculation and 
(|23| ) stays the same with these replacements. The "classical" mass terms T B (k = 0) 
of the bosons B are given by Att 2 . 

4.2 Loop corrections to the bosonic propagators 

We can now insert the explicit formulae for the vertices from appendix |C| and eval- 
uate the inverse propagators for the various boson species. The boson p corresponds 
to the charge density p ~ ipl^a (see appendix |B| for further details). In the basis 
(Pi,p 2 ,P3>Ai) one finds the correction 

AT^(K) = AhlT r ^ 



diag{ [cos(fc 1 /4) cos(A; 2 /4)] 2 {<£(+, +) + <£(-, -)} , 

- [sin^/4) sin(k 2 /A)} 2 {g r + (+, +) + g r + (-, -)} 
[oos(*i/4) sin(fc 2 /4)] 2 {<£(+, -) + g r _(-, +)} , 

- [sin(^/4) cos(A; 2 /4)] 2 {<£(+, -) + g r + (-, +)} 

Here we have defined the Matsubara sums (m, n G Z, e,; G {+, — }): 

r,e ( \ _ \^ (ci + eic 2 )(c / 1 + e 2 c 2 ) + e 3 uu' 
9e 3 tei, e 2 j - ^ [( ^ + ^ + ^ [(c , + ^, )2 + ^ 

, Si(m,a ei ,& e2 ) S 2 {m,a ei ,b e2 ) 
= (ci + e lC2 )( Cl + e 2 c 2 ) ^ ± e 3 ^ 



(26) 



(27) 



which depend on m, k and g. The upper sign applies to real bosons (r) and the lower 
sign for complex bosons (c). The sums Si j2 can be found in eq. (|50| ) in the appendix 
|A[ The frequencies u, uo' appearing in the definition of g read 

lu = 2ir(n + 1/2)T ^' = / 27r ( m + ra + |) T for real bosons . . 

\ 2n(m — n — |)T for complex bosons 

For the arguments of Si j2 we use the abbreviations (e^ G {1,-1}) 

a ei = (ci + eiC 2 )/(7rT) 

6 £i = (c' 1 + e iC2 )/(7rT), (29) 

where q and c[ are given by 

/ /n \ i \ 2tcos((ki + Oj)/2) for real bosons ,„-,. 

q = 2tcos(qi2). d = < _. )), ^V'.' , , . . (30) 

[ 2r cos((kj — qi)/2) tor complex bosons 
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From eq. fl26[ ) it is obvious that only p 1 receives a propagator correction for k = 0. 
This reflects the fact that only p l has a nonvanishing vertex in this limit. 

A similar formula is found for the boson p, if one replaces 

cos(fcj/4) <-> sin(fcj/4), /i p — > /i p . 
Similarly, the expressions for q x are obtained by 

cos(fti/4) <-> sin(fei/4), h p — > h qx 
and for q y by 

cos(fc 2 /4) <-> sin(fc 2 /4), h p -» /i 9w . 

In the case of the bosons with spin index, m,a,~g xy , one obtains the same result as 
for the above scalar bosons p, p, q x , q y since tr (tjTj) = 25^ = 5ytrl 2 spm . 

The bosons s, c, £ again receive a loop correction with the same structure as the 
corresponding scalar ones because the vertices are the same - one only has to mul- 
tiply by —2 and to replace c[, c 2 and u' by the corresponding expressions for the 
complex bosons (J30|) , (p7j) . 

For the boson d one finds: 

/2, 



d d J-« (2tt) 2 ^ 



diag{[cos(A;i/4)cos((A ;2 - 2g 2 )/4) - ecos(A; 2 /4) cos((fci - 2 gi )/4)] 2 ^(e, e), 

- [sin(A;i/4) sin((A; 2 - 2g 2 )/4) - esin(A; 2 /4) sin((fc 1 - 2 gi )/4)] 2 ^(e, e), 
[cos(A;i/4) sin((A; 2 - 2g 2 )/4) - esin(fc 2 /4) cos((fci - 2q l )/A)] 2 g c L{e 1 -e), 

- [sin(A;i/4) cos((A; 2 - 2g 2 )/4) - ecos(fc 2 /4) sin((A; 1 - 2 gi )/4)] 2 ^_(e, -e 

(31) 
From this we obtain the inverse propagator for the boson e by the replacements 
g e3 (e 1 ,e 2 )-*g e3 (-e 1 ,-e2) } h d ^h e 

Similarly, the propagator correction for the boson v x is 

AT^\K) = ~2hl T r -j% 
Vx Vx J_ n (2vr) 2 

diag{[sin(fci/4) cos((A; 2 - 2g 2 )/4)] 2 {g c _(+, +) + «£(-, -)} , 



- [008(^/4) sm((k 2 - 2g 2 )/4)] 2 {<£(+, +) + g c + (-, -)} 
[011(^/4) sin((A; 2 - 2g 2 )/4)] 2 {<£(+, -) + g c _(- +)} , 

- [008(^/4) cos((A; 2 - 2g 2 )/4)] 2 {<£(+, -) + <£(-, +)} 



(32) 
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and for the boson v y one obtains 



^{K) = -2hlT 



d 2 q 



-* (2vr) 2 
diag{[sin(A; 2 /4) cos((h - 2 gi )/4)] 2 {<£(+, +) + «£(-, -)} , 

- [cos(A; 2 /4) ^((fcx - 2 gi )/4)] 2 {£(+, +) + «£(-, -) } , 
[cos(A; 2 /4) cos((fci - 2 gi )/4)] 2 {«£(+, -) + g c _(-, +)} , 

- [sin(A; 2 /4) sin^ - 2 gi )/4)] 2 {<£(+, -) + <£(-, +)} j. 



(33) 



Though it is not always immediately apparent, the above expressions for 
Ar^ (X, X') are symmetric under reflection of the external momenta. Note fur- 
thermore that the sums Si are symmetric in the Matsubara frequency, Si(m, a, b) = 
Si(—m, a, b). 

Up to now, we did not consider the off diagonal propagator terms involving bosons 
of different species. Many of these terms are zero due to symmetry arguments: The 
real and the complex bosons do not mix because of the t/(l)-symmetry (which can 
be interpreted as charge conservation) and the charge waves do not couple to the 
spin waves because of the SU (2)-sjrametrj acting on spins. Unfortunately, there is 
no symmetry prohibiting the mixing of the site and link pairs. The best thing one 
can do is considering the transformation 



(x,-r) (34) 



which is no symmetry of the original fermionic partition function for p 7^ (since 
this transformation yields a theory with a chemical potential with reversed sign). 
Nevertheless it is a symmetry of the fermionic effective action, since the latter does 
not depend on the chemical potential (this transformation reverses the sign of p and p 
and correspondingly the sign of the associated source L' -\-p). To keep this symmetry 
in the partially bosonized effective action, we find the corresponding transformation 
properties of the bosons 



(fr\ 




f ^ ^ 




(w\ 




( -M \ 


^2 
^3 


(x , r) -► 


-^2 
^3 


(x,-t), 


r 2 
r 3 


(x, r) -> 


r 2 


Ww 




\-^J 




\rj 




\ n J 



$(f,r) 



$(x, 



-T 



X(x, T ) 



■X(x, 



-T 



Formulating this in Fourier space, we find that this symmetry tells us that the prop- 
agator matrix elements for mixing between site and link pairs are uneven functions 
of the Matsubara frequency m. In particular, for m = these propagator matrix 
elements vanish. 
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Furthermore, we find by explicit calculation that different colors do not mix even 
for different boson species. The only remaining off diagonal terms in the inverse 
propagator occur between charge waves with equal colors and correspondingly spin 
waves with equal color as well as site pairs and link pairs with equal color. It is 
tempting to diagonalize further in these remaining non diagonal 4x4 propagator 
blocks. For nonvanishing momentum this is a highly nontrivial task und subject to 
current work. Suffice to say that nearly all the non vanishing offdiagonal terms which 
couple different species vanish in the limit of homogeneous fields. To be precise, only 
the link pairs of different species couple even in the limit of homogeneous fields. 

4.3 Mean field results 

In H we have evaluated the fermionic loop correction to the effective action for 
fields constant in time and space, thus obtaining the effective potential. In par- 
ticular we have analyzed the bosons p 1; H2 and d\. If one is mainly interested in 
antiferromagnetic and superconducting behavior we have now shown a posteriori 
that it was legitimate to choose exactly these color-combinations. They are stable 
under fluctuations in the sense that the minimum of r^ 2 ^ is at K = and they have 
a nonvanishing coupling to the fermion bilinears. 

To make the connection between the present results for the bosonic propagator and 
the mean field calculation of more explicit it is interesting to investigate the 
above expressions for constant fields, i.e. in the limit of vanishing outer momenta. 
In this limit, we should have T ( a\K = 0) = M 2 = 2 |^| =d=g= =Q , where U is 

the mean field potential defined in [0 and correspondingly T d (K = 0) = Mj = 

8Uo \ 

dd? I p=d=a=fi=0' 

Indeed, as already anticipated in figure |], the results are equal and we find 

Ml = * - 2h l f *L E ^ (*(*+"»>) (35) 

2 2 1 2 [" d2( l ST tanh(^(ci + ec 2 )) 2 

* = 4 * -l^LwJ^, ^T7^ < c '-« 2) ' 

In a similar way, mass terms for the whole set (cf. appendix |B|) of bosons have been 
calculated in this limit. In fig. (|5]) we show the mass of different bosons as a function 
of temperature. We choose h 2 = h? p = h\ = h 2 , (the couplings for the other bosons 
are then uniquely determined by the conditions in appendix |FJ). Only those bosons 
are shown which reach zero mass at T/U > 0.01. It is reassuring to see that the 
most prominent degrees of freedom as ferromagnetic, antiferromagnetic, s- and d- 
wave superconducting states, which are known to play a role in the Hubbard model 
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near half filling, emerge naturally in our framework as the candidates competing in 
determining the way of symmetry breaking. On the other hand, we see that because 
of the arbitrariness of the Yukawa couplings we cannot decisively identify the boson 
that wins the mass run to zero. 

Getting rid of the arbitrariness of the Yukawa couplings will be a significant premise 
for correct predictions of phase transitions. This is also necessary for a quantita- 
tive comparison with other methods. Qualitatively, our results are compatible with 
renormalization group investigations of the four fermion interactions in the Hubbard 
model [|J. We stress that our formalism for partial bosonization is exact and the final 
result has to be independent of the choice of the Yukawa couplings. This requires the 
inclusion of the bosonic fluctuations which we have omitted in the present work. Also 
note that in our figures the bosons are treated as completely independent. The mass 
of every boson is calculated for a vanishing expectation value of all bosonic fields. 
Going beyond this assumption, one has to take into account the fact that in case of a 
transition to a phase with a nonvanishing order parameter, the corresponding boson 
no longer will have a zero expectation value. Often this nonzero expectation value 
will prevent the masses of the other bosons going to zero — a behavior encountered 
in [0 for the competition between antiferromagnetic and superconducting phases. 



5 The high temperature limit 

One of the goals of the one loop calculation was to find a useful ansatz for the bosonic 
propagator that can be used in a renormalization group study of the Hubbard model. 
For T — > oo, the fermion loop corrections vanish and we end up with the classical 
bosonic mass terms. Letting T — > 0, we have to face the complications of divergencies 
at the Fermi surface. Consequently, we have to lower T from oo to some temperature, 
that is low enough to reveal nontrivial physical behavior and high enough to keep 
this behavior so simple that we are able to write down well justified, easy analytical 
approximations to the one loop expressions at this temperature. This aim in mind, 
we calculate a high temperature expansion of our one loop expressions, keeping as 
many terms in an expansion in T _1 as is necessary for the result to be nontrivial. 

Because in our expressions the sum 5*1 always occurs with a factor T~ 4 and S 2 with 
a factor T -2 , it will suffice to expand S 2 up to order T° if we want the result to be 
correct up to order T™ 2 , which will turn out to be the lowest interesting order. We 
find to this order only a contribution for m = (^(m 7^ 0, a, b) = 0(T~ 2 )) 

- 2 -.-_„, 



resulting in 



S 2 (0,a,b) = - + O(T- 2 ) (36) 



^(ei,e 2 ) - s 2 T) 2 ' 0«s( Cl ' e2 ) ~~ ~ (2T) 2 ' ^ 
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Figure 5: Mass terms of different bosons as a function of temperature. The upper 
plot is for h 2 /U = 10, the lower plot for h 2 /U = 40. 
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We will mainly be interested in bosons that have a nonvanishing loop-correction in 
the limit k — > and thus only list those below. In the high temperature limit one 
thus obtains: 

Alf^) = __^[cos(A;i/4) cos(A; 2 /4)] 2 cW (38) 

The same result applies — with appropiate replacements of the Yukawa couplings 
- to p 2 ,Q x ,4: an d 9y,3 and hence also to mi,a2,g X 4,g y3 . Noting the extra minus sign 
in g c relative to g r in the high temperature limit this also, apart from a factor of 
+2, applies to the result for the complex bosons si, c 2 , t x ^t y ^ • 

For d and e one finds 



Ar SW = ~^m0 



h 2 

diag{[cos 2 (A;i/4) + cos 2 (A; 2 /4)], [sin 2 (fci/4) + sin 2 (A; 2 /4)], (39) 

[cos 2 (A;i/4) + sin 2 (fc 2 /4)], [sin 2 (A;i/4) + cos 2 (A; 2 /4)] 



Finally, for v x ^ v x> ± v y $ and v y ^ we find solutions with "stripes" in the x- or y- 
direction 

h 2 
,(2) ri. 






(40) 



Wj/,23^ I 2T 



We want to stress that nearly all expressions which give a contribution in the limit 
k — > (except 62,3,4, ^2,3,4) possess minima at k — 0. This means that all physically 
important bosonic degrees of freedom have the tendency to prefer homogeneous field 
configurations. 

As only the term in the Matsubara sum with m = receives a correction in the high 
temperature limit, we define 

Ar (2) (fc) = Ar (2) (cu m = 0,£). (41) 



The high temperature limit for the a and d bosons coincides with good accuracy 
with the graphs shown in fig. (0). 
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6 Exact renormalization group 

The inclusion of the bosonic fluctuations and the exploration of small T can pre- 
sumably best be achieved by using non-perturbative flow or renormalization group 
equations for the scale dependence of the effective action. This short section ar- 
gues that the results of the present work may constitute a good starting point both 
for the formulation of the flow equations and the setting of "initial conditions" at 
microscopic scales. 

For the formulation of an exact renormalization group equation in the framework of 
the effective average action || we add to the action (§) an infrared cutoff 

A k s = J2{^*(Q) R kF(Q)i>(Q) (42) 

Q 

+ J2%(Q)R7b(Q)%(Q) 4E^H)€(*W) 

7 

where we adopt a vector notation for the fields, making the color indices implicit. 
Accordingly, R k F(Q) is a matrix in spin and color space, and R k B(Q) is a matrix 
in color space. The effective action T transmutes to the effective average action T^ 
with r fc=0 = T if R k=0 = 0. 

For the bosons we suppress the low momentum modes byQ 

R kB {Q) = Z B , k X k (Q)9(X k (Q)) (43) 

Here the form of X k (Q) can be adopted to the kinetic term in the inverse bosonic 
propagator 1^ such that T B ' k + R Bk = M k + k 2 becomes independent of Q for 
small Q 2 . The precise form of X k (Q) and the wave function renormalization Z Bjk 
generically differ for the various bosons. 

For the fermions we want to "regularize" the Fermi surface and therefore use an 
infrared cutoff that guarantees that the square of the inverse average propagator 
is always positive and different from zero. The cutoff then acts like a gap in the 
fermion spectrum. Since the temperature has the desired properties, we propose 

R kF (Q) = 2mnkZ F>k . (44) 

This cutoff simply replaces the temperature T in the inverse fermionic average prop- 
agator T —>■ T + k. 

The exact evolution equation for the scale dependence of T k reads 

d k T k = ^STr{d k R k (r^ + R k )-i} (45) 



10 This cutoff resembles the opimized cutoff |J. 
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where r^ 2 ^ is the matrix of second functional derivatives with respect to both bosonic 
and fermionic fields. The supertrace STr involves a momentum integration as well 
as sums over Matsubara frequencies and internal indices. It contains a minus sign 
for the fermions. The cutoff Rk is a block diagonal matrix, where each block is 
equal to one of the infrared cutoff functions introduced above. The inverse average 
propagator (T k + Rk) becomes particularly simple for our proposal (|3|), (fH|). 

Due to the simple form of the fermionic cutoff we may use the results of this paper 
for a one loop calculation of Tk for high k. This simply replaces the prefactor T _1 in 
AI^ 2 -* by T/(T + k) 2 in the results of sect. 5. For large enough k the approximation 
should be reliable and we can use this result as an "initial value" of the (functional-) 
differential equation P5| ) for large k. The aim will then be an approximate solution 
for k — > whereby the effective action is recovered. 



7 Discussion and conclusions 



The colored Hubbard model [|| is an equivalent reformulation of the usual Hubbard 
model in two dimensions as a Yukawa theory. The bosonic degrees of freedom cou- 
ple to fermion bilinears which are chosen such that they reflect the most interesting 
physical properties of the system, i.e. antiferromagnetic or superconducting behav- 
ior. This has the advantage of representing these properties explicitly as bosonic 
expectation values instead of dealing with complicated properties of fermionic inter- 
actions. 

A previous mean field analysis suggests the occurrence of spontaneous symmetry 
breaking in the antiferromagnetic or d-wave superconducting channels for low T. The 
present computation supports a crucial ingredient for such a mean field calculation, 
namely that the composite bosonic condensates are spatially homogeneous. 

In this article, we have calculated the one loop corrections to the bosonic propagators 
analytically. The remaining momentum integrals cause trouble for low temperature, 
where the singularities near the Fermi surface come into play. On the other hand, 
we are able to perform the momentum integrals in a high temperature limit. We 
find that the kinetic properties of the bosons emerge in order T~ 2 . In this order, 
the momentum dependence of the propagators is extremely simple and shows a 
number of pleasant properties. The corrections to the inverse propagator T^ are 
negative definite, which means that they do not change sign when the parameters 

(2) 

are varied. Furthermore T B becomes minimal for homogeneous field configurations 
for all interesting boson species, which means that these field configurations should 
be preferred by the system. This result yields a strong argument in favor of the 
basic assumption of the mean field approximation in 0. Our task to extract the 
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most interesting degrees of freedom of the Hubbard model from the complicated 
momentum structure of its vertices has been accomplished to some extent. 

A more accurate treatment of the low temperature properties of the colored Hub- 
bard model may use exact renormalization group equations to analyze the flow of 
the Yukawa couplings and the effective potential. This should reveal the occurrence 
of spontaneous symmetry breaking in certain ranges of the density and the temper- 
ature ||. A straightforward approach to this task is complicated by the problem to 
write down a suitable truncation for the effective average action, which on one hand 
should be simple enough to be tractable, but on the other hand should contain the 
interesting physical behavior. A common way to compromise these two demands 
is a truncation which copies the classical action, merely introducing wave function 
renormalization constants, flow dependent couplings and some terms for the effec- 
tive potential in agreement with the symmetries of the system. Unfortunately, for 
the colored Hubbard model this approach is too simple, because the bosonic prop- 
agators of the classical action consist of mass terms only. Some kinetic behavior of 
the bosons obviously has to be included in the truncation. A priori it is far from 
clear in which way this should be done. We pursue here the way to calculate the one 
loop corrections of the bosonic kinetic terms to get an impression of how the bosons 
behave once fluctuations are included. 

We finally want to comment briefly on the formal status of our computation as 
compared to the first order in a perturbative expansion of the Hubbard model in 
a purely fermionic language. On the one hand our computation goes beyond the 
perturbative result for the effective four fermion coupling. This is apparent from 
eq. fll8|) if we realize V ~ U 1 ^ 2 and T B = M + cU. An expansion of T B yields a 
contribution to the ^-interaction ~ U 2 similar to first order perturbation theory 
but also higher order terms which amount to a resummation. In particular we find a 
divergence of the four fermion interaction for a critical temperature (when T B van- 
ishes) and the onset of instability towards spontaneous symmetry breaking. These 
features cannot be seen in perturbation theory. On the other hand, our calculation 
does not reproduce the complete perturbative result ~ U 2 . This would require in 
addition the computation of V in order U 3 ^ 2 as well as of four fermion interactions 
~ U 2 which remain lPI-irreducible even in the partially bosonized language. Our 
computation of the effective four fermion vertex can therefore be trusted only if the 
selected channels really dominate as suggested by the low temperature behavior. 

The present insufficiency of our setting is also apparent from another perspective. 
Our results in order U 2 depend strongly on the choice of the Yukawa couplings h 
whereas no such parameter appears in the perturbative calculation in the fermionic 
language. This severely limits the predictive power of our appoach so far (cf. fig. ^, 
[5]). The present one loop calculation may be used for a determination of h by an 
optimization procedure. In fact, one may compute the four fermion vertex in selected 
channels (e.g. antiferromagnetic, d-wave-superconducting and charge density) both 
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by our "one loop result" and by perturbation theory. Since the "one loop result" 
depends on three Yukawa couplings those may be fixed by matching the term ~ U 2 
to the perturbative result in these channels. This guarantees that the neglected 
corrections to h and lPI-vertices are indeed small for the selected channels. 

In consequence a mean field or renormalization group calculation with such an op- 
timized choice of h would start close to perturbation theory for high T (at least for 
the selected channels) and nevertheless go far beyond at low T. One may then hope 
that further corrections in the other channels may only result in minor modifications. 
We hope that a combination of such an optimization of the choice of Yukawa cou- 
plings together with the study of non perturbative renormalization group equations 
may permit a reliable quantitative computation of the properties of the (colored) 
Hubbard model. 



Appendix 



A Useful Formulae 
A.l Abbreviations 



Q = (uJ n ,q), X~(t,x), 

Z bosons 



QX = uo n T + xq, uo n = 2imT, n G -, „ .. , ,, 

' ' Z + 1/2 fermions 

X J 3 Q n J-^ Z7[ ) 

S(Q - Q') = ^6 n>n , ■ (2n) 2 5(q- q') 
S(X - X') = 5{t - r') • S(x - x') (46) 

Note that S(q — q') is periodic in 2n i.e. S(q + 2i\e^) = S(q). The same applies to 
S(t) = ±<5(r + (3) for bosons/fermions. 
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A. 2 Fourier transforms 

For the transforms of the fermionic if), ip* and bosonic fields B, B* (generically 
denoted by x, ■%*) we use 



(47) 



All 



* a (X) = v^^exp (-i{QX + z a q}) xKQ) 



a a\ _ fa a 



l ^ l ~2'2>'' ^ 2= V2'2/ 

a a\ ... / a en ^ 



^ _ U'~2>/' * 4 ~V 2' 2 
We set the lattice distance of the coarse lattice (c.f. fig. [l|) to unity in our calculations: 
2a = 1. 

As a consequence of the choice of z a we find simple transformation properties of 
the fermionic Fourier modes with respect to the translations by a in the x- and 
^-direction, T x and T y . With s x = (o,0), s^ = (0, a), ip(Q) = (i/>i(Q),... ,i}>i(Q)) 
one obtains 

TMQ)=^( T i °)v>(Q) 

^(Q)=c*(° q)^(Q) (49) 



and similarly for ip*(Q) with e l?5 replaced by e 



iqs 



A. 3 Some sums 

The following sums are useful when evaluating the Matsubara sum appearing in the 
one loop corrections to the bosonic propagator: (uj m = 2irmT, m G Z) 

Si (m, a, b) := 2_) 



z [(2n + l) 2 + a 2 ] [(2(n + m) + l) 2 + 6 2 ] 



7r 6(4m 2 - a 2 + 6 2 ) tanh(^) + a(4m 2 + a 2 - 6 2 ) tanh( 



7rbA 
2 I 



S 2 (m,a,b) := ^ 



ab[Am? + (a + b) 2 ] [Am 2 + (a - b) 2 } 

(2n+l)(2(n + m) + l) 



■ 



(50) 



new, 



[(2n + l) 2 + a 2 ] [(2(n + m) + l) 2 + 6 2 ] 



vr a(4m 2 + a 2 - 6 2 ) tanh(f ) + 6(4m 2 - a 2 + 6 2 ) tanh(f 
2 [4m 2 + (a + 6) 2 ][4m 2 + (a - 6) 2 ] 
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(51) 



B Fermion bilinears and bosons 



B.l Naming scheme 



Neutral real bosons 


charge waves R 
spin waves S 


P, P, Qx, % 
™>, a, g x , g y 


Charged complex bosons 


site pairs $ 
link pairs x 


$j C, Z XJ by 
6, Utj V X i Vy 



All these boson appear in four distinct colors. The corresponding fermion bilinears 
are designed by a tilde, e.g. p, p, etc. 

B.2 Fermion bilinears 

For each of the bosons b(X) introduced above, we have a corresponding bilinear 
b(X) which couples to the boson b(X). Suppressing the dependence on X, we define 



0ab 

Xab 

~* 
Xab 



Vb^a 

i)l{ir 2 )^a 
-ipZ(iT 2 )ip* 



(52) 



and the composite bilinears 



P = #11 + #22 + #33 + #44 

P = #11 — #22 + #33 - #44 

% = #11 + #22 - (733 - #44 

<ix = CTll - #22 - 0-33 + #44 



m 

a 

9y 

9 X 



Hi + 022 + 033 + 044 
>11 - 022 + 033 - 044 
hi + 022 - 033 - 044 
hi ~ 022 - 033 + 044 



S = Xll+ X22 + X33 + X44 
C = Xll - X22 + X33 - X.U 



t 



!l 



Xll + X22 - X33 - X44 
Xll - X22 - X33 + X.U 



e = Xl2 + X23 + X34 + X41 
d = Xl2~ X23 + X34 - X41 



V 



>l 



Xl2 - X34 

X23 - Xn 



(53) 



Note that q x / y , g x / y , t x / y are linear combinations of the definitions given in 0. 
With these new definitions, all bosons have a simple transformation behavior under 
translations by a. 
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B.3 Yukawa couplings 

Integrating out the bosons in ([?[) gives a purely fermionic theory. However, this 
theory coincides with the Hubbard model only under certain conditions for the 
Yukawa couplings. These conditions are (with h% = ^H^U) 

H p = 3(A 2 - A 3 ) H* = 2Ai + A 2 + 3A 3 + 1 

H p = 3(A 2 + A 3 ) H s = 2Ai + A 2 - 3A 3 + 1 (54) 

H q x = H q y = 3A 2 H Sx = H Sy = 2Ai + A 2 + 1 



3 

H s = H c = H tx = H ty = -Ai, 2H e = 2H d = H Vx = H Vy = 6A3. (55) 



where the parameter Aj obey 

A 4 >0, i = 1,2,3, 

A 2 > A 3 , (56) 

2Ai + A 2 + 1 > 3A 3 . 

B.4 Bosonic colors 

The color index for the fermion bilinears and the bosons is defined by 

w ly (x) = T y T~ 1 u) 7 (x) , W2 1 (x) = TyW^(x) , 

W3 7 (x) = W 7 (x) , W4r/(x) = T~ 1 W.y(x) (57) 

and similar for u,w,u, where the fermion bilinears without color index are given by 
the definitions in sec. (|B.2|) . 



C Vertex factors for the Hubbard model 



The vertices V W (Q',Q") for the bosons p, p, q x>y depend only on the momentum 
Q = Q' — Q". With e x = (1, 0), e y — (0, 1) and z a , a — 1 . . . 4, given in the appendix 
0, eq. d48|) , they can be written in the form 

V£ C (Q', Q") = V^ C (Q) = ^e-^e^M^Q) ® lf n , (58) 
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The color matrices M™ read 

M?{Q) = diag{l, e iS *\ e^*'^, e'^}, M P (Q) = diag{l, 1, e~ iS ^, e'^}, 

M§(Q) = diag{l, 1, 1, 1}, M^Q) = diag{l, e** e** 1}; 

M?(Q) = (-l)- 1 diag(l, -1, 1, -1) M?(Q); 

M! y (Q) = Mf (Q) • diag(-l, -1, 1, 1), M*(Q) = M£(Q) • diag(-l, -1, 1, 1), 

M^(Q) = Mg(Q) ■ diag(l, 1, -1, -1), Ml y {Q) = M£(Q) ■ diag(l, 1, -1, -1); 

Mf*(Q) = Mf (Q) • diag(-l, 1, 1, -1), M%*(Q) = M P 2 {Q) • diag(l, -1, -1, 1), 

Mt(Q) = Mg(Q) • diag(l, -1, -1, 1), M|»(Q) = M£(Q) • diag(-l, 1, 1, -1). 

(59) 

The same can be obtained for the bosons with spin index, rh, a, g XjV , by substituting 

j spin _^ z^pin 

For the bosons s,c,t X}V one finds similarly (c spin = ir<i): 

V£ C (Q',Q") = ^e^W%^W')M£ tC (Q',Q'')®J** a , 



Mf(Q',Q") = M?(-Q' - Q"), M?(Q',Q") = M?(-Q' - Q"), 
M$ (Q\ Q") = Mf (-Q' - Q"), Mf (Q', Q") = Mf {-Q' - Q"), 



(60) 



while d, e, v x>y are a bit more complicated. Let us define e JJ = e^ Ziq +Zjq •* and a 
♦-product C = A*B by CV,- := AijBij (no sum over indices here!). One then obtains 

Vf(Q',Q") = !^e- ig °® +ir >M?(Q',Q")®d* in 



( o 
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r -'( "' ' « 
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e e 







e 14gi(j e a 
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e 32 e i,fe H q e 34 e 

e 41 e i<fe H g e ^ e i(q , +q")e y Q 
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li p iq"e y \ 
/ 



„21 -iq'Sx 
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\ e 4i e 43 ) 









e 12 e -iq"e 



e 32 e -i(q'+q")e. 
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e 34g-ig'e 



e iZp-iq"e x 
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(61) 
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With the aid of the ^-product the other vertices can now be obtained from these 

/ 1 -1 \ 

Mf(Q',Q") = 



M V C *(Q',Q 



M v c y (Q' } Q") 



10-10 

0-101 

\ -1 1 / 

/ -1 \ 

10 

10 

\ -1 / 

/ 1 \ 



* Ml (Q', Q' 



M c e (Q',Q")-A c , A = (-1,1, 1,-1); 



10 

0-1 

\ -1 / 



* Ml (Q' } Q") • A c , A = (-1,-1, 1,1). (62) 



The transition from the ■u*^'0-vertices V u * to the w^^-vertices V u can be carried 
out by 

V U *(Q',Q") -> V^(Q',Q") = -(V u *(g',Q"))* = -V u \-Q\-Q"). 
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